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Abstract: 

We discuss dark matter decay mediated by exotically charged particles ("exotics") in a 
supersymmetric model with two dark matter (DM) components: One is the (bino-like) 
lightest supersymmetric particle (LSP) Xi and the other is a newly introduced meta-stable 
neutral singlet A^. decays to x^'^^~ via a dimension 6 operator induced by a penguin- 
type one loop diagram with the life time of 10^^ sec, explaining energetic cosmic excess 
observed recently by PAMELA and ATIC/PPB-BETS. The superheavy masses of exotics 
(~ 10^^~^^ GeV) are responsible for the longevity of A^. The superpartner of develops 
the vacuum expectation value (VEV) of order TeV so that the DM A^ achieves the desired 
mass of 2 TeV. By the VEV, the U(1)r, symmetry is broken to the discrete Z2 symmetry, 
which is identified with the matter parity in the minimal supersymmetric standard model 
(MSSM). Since we have the two DM components, even extremely small amount of A^ 
[0(10^"'^'^) < {riN/n^)] could account for the observed positron flux with relatively light 
exotics' masses [lO^^ Q^y < j^f^^^ < iqI^ GeV]. 
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1. Introduction 

The recently reported observations by PAMELA ||, § and ATIC/PPB-BETS g, g col- 
laborations on excess of high energy positrons from cosmic ray have attracted more and 
more attentions. As many literatures pointed out, dark matter (DM) decay H, |6|, or 
annihilation |§ would be deeply involved in the observed positron excess.^ If it is indeed 
caused by DM, however, the observation should be accepted as a puzzle, because it is hard 
to be understood within the framework of the conventional DM scenario, particularly, by 
the minimal supersymmetric standard model (MSSM). Thus, the observed positron excess 
might be a hint toward a new physics beyond the standard model (SM). 

As noticed in Refs. |^ , the positron flux needed to explain the observation of ATIC /PPB- 
BETS (and also PAMELA) can be produced by leptonic decay of DM [|, with 2 
TeV mass (= w-dm) via a dimension 6 operator (four fermion interaction) suppressed 
by MQyrp ~ (10-*^^ GeV)^, by which the decay rate is estimated as 
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Hadronic decay channels should not exceed 10 % to be consistent with the PAMELA'S 



data Q. The DM decay scenario avoids the constraint from the 7 ray flux [10| by the 
HESS observations of galactic ridge [|l^]. However, it is not trivial to see which physics at 
the Mqut scale can provide such a low energy effective four fermion interaction, allowing 
DM to decay dominantly into the SM leptons: In most of grand unified theories (GUTs) 
embedding the SM, the gauge interactions by superheavy gauge boson exchanges can easily 
provide four fermion interactions suppressed by Mq^^^. But they do not prefer only such a 
leptophilic decay mode of an electrically neutral particle or DM. Hence, one should explore 
the possibility of a leptophilic Yukawa interaction for DM decay. 

Recently, supersymmetric (SUSY) models possessing one more dark matter component 
N apart from the (bino-like) lightest supersymmetric particle (LSP) x have been suggested 
as the resolutions of the PAMELA/ ATIC anomaly ||l^, ^. Particularly in the model 



""^Alternatively, astrophysical sources such as pulsars could explain the positron excess 
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of Ref. 1^], the anomaly is explained by decay of the extra DM component A'^ into x^^e^ 
through a dimension 6 operator. The effective dimension 6 operator for DM decay is ob- 
tained from some renormalizable leptophilic Yukawa interactions with the dimensionless 
coupling of order unity, after a pair of vector-like SU(2) lepton doublets {L, L'^) and lepton 
singlets {E, E'^) decoupled. The superheavy masses of lW, (~ IQi^ GeV) are respon- 
sible for the longevity of N. Since the gauge group is just that of the SM and the low 
energy field spectrum is the same as that of the MSSM except the neutral singlet N, the 
gauge coupling unification in the MSSM is protected in the model. This model is easily 
embedded in flipped SU(5), which is a leptophilic unified theory |13|. 

Most of phenomenologically promising string models predict a lot of vector-like super- 
heavy exotic states ("exotics") carrying fractional electric charges [|l4|. One might expect 
that such superheavy exotics also can play the role of (L, L'^) and (E, E^) in the model 
of Ref. 0, mediating DM decay via the dimension 6 process. Their superheavy masses 
(~ 10^^ GeV) could lead successfully to lO^^ sec. life time of the DM as desired. Con- 
sidering the case that superheavy exotics mediate DM decay, however, one should notice 
a remarkable point: Most of all, fractionally charged heavy particles can not decay to the 
light SM leptons, because of the charge conservation. Thus, if exotics are involved in the 
process, ^ e"*" -|- e~ -|- neutral particles, where the initial and final states are the states 
only with the integral electric charges, they should be co-created and co-annihilated be- 
tween the initial and final states. It means that DM decay is possible only at loop levels, 
if exotics dominantly mediate DM decay. 

In this paper, we explore the possibility that DM decay is mediated by a one loop 
diagram. If the mediators are indeed fractionally charged superheavy particles, we should 
necessarily consider the loop induced process. However, our study is not confined only to 
the case of fractionally charged heavy field mediation, but covers more general cases of 
loop induced DM decays. 



2. The model 



Let us consider the vector-like superheavy superfields {E,E^), (X,X^), and {0,0^). Their 
quantum numbers are shown in Table 1. If g is a fractional number, E^'^\ X^'^\ and O^^^ 
become regarded as exotics. In Table 1, we present only the first generation of the charged 
lepton singlets, e^. Concerning the R charges of the other MSSM superfields, we assign 
1 to the MSSM matter superfields hke and to the two MSSM Higgs doublets. We 
leave open the possibility that E^'^\ X^'^\ and O^^^ are charged also under other (visible 
or hidden gauge) symmetry Q. For the case that this model is embedded in flipped SU(5) 
[=SU(5)xU(l)x], G can correspond to SU(5). 

If the deviation of (e"*" -|- e~ ) observed by ATIC /PPB-BETS from cosmic ray is indeed 
caused by DM decay, the mass of DM should be around 2 TeV Q. In order to protect 
the status of SUSY as the solution of the gauge hierarchy problem, we should assume that 
the mass of the LSP is of O(IOO) GeV or lighter. Apart from the (bino-like) LSP X) thus, 
we introduce one more dark matter component with 2 TeV mass, which is the fermionic 
component of N in Table 1, to account for the ATIC/PPB-BETS' data. 
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Table 1: The hypercharges and R charges of the superfields. The hypercharge q can be a fractional 
number. The vector-like exotic superfields, E^'^\ X^'^\ and O'^'^^ are all decoupled from low energy 
physics due to their heavy masses. The (visible or hidden) symmetry Q is optional. 

The relevant superpotential in our model is composed of the trilinear and bilinear 
terms: W = Wtri + Wbi, where Wtri and W\^\ are, respectively, given by 

Wtri = NEX + XOe" + iV^ (2.1) 
Wu = MeEE" + MxXX^ + MoOO^. (2.2) 

We dropped the dimensionless Yukawa coupling constants in Eq. ( p.lD for simplicity. They 
are tacitly assumed to be of order unity. The dimensionful parameters. Me, Mx, and Mq 
in Eq. {^) are lO^^-lO^^ GeV. Thus, the vector-like fields {E,E''), (X,X'=), and [0,0") 
are superheavy. To avoid couplings with the other charged lepton singlets, and t^, 
one can introduce a family dependent U(1)pq symmetry. It can explain the smallness of 
the electron mass m, g. The term in Eq. ( ^.l| ) is introduced such that the scalar 
component of N, i.e. N promptly decays into the two fermionic components 2N: The 
mass of the fermionic component of A'^ (~ 2 TeV) is induced by the vacuum expectation 
value (VEV) {N). On the other hand, the mass squared of the scalar component of is 
given by (|(iV)p + "^-3/2)5 where rn'^^^ comes from the soft scalar mass term of N. We will 
discuss later how the VEV of N could be developed. We just assume that the soft mass of 
N is heavy enough (> 4 TeV) for the decay 2N to be possible. Since we don't want 
the N'^ term with a too large mass parameter in the superpotential, we employ the U(1)r 
symmetry to forbid it from the bare superpotential. 

This model is easily embedded in flipped SU(5) [^]. To account for the PAMELA'S 
important observation, i.e. no excess of anti-proton the lepton singlet should not 
be accompanied with quarks in Eq. ( |2.1| ), when the model embedded in a GUT. Since 
in flipped SU(5) e'^ and A^ remain SU(5) singlets, I5 and lo, respectively, flipped SU(5) 
models can be perfectly consistent with the PAMELA'S data [|l^. Moreover, flipped SU(5) 
is phenomenologically attractive: The notorious doublet/triplet splitting problem in GUTs 
is very easily resolved via the missing partner mechanism |]l5| . The predicted fermion mass 
relation in flipped SU(5) is just that between up-type quarks and Dirac neutrinos masses. 
Since the Majorana neutrino masses are still not constrained, however, the mass relation 
in flipped SU(5) does not encounter any difficulty in matching the real data on fermion 
masses. 

The presence of the A-term corresponding to A^^, i.e. {my2^'^ + h.c), and |iV|^, 
(and also the soft mass term 771^^2 |A^P) the scalar potential permits two vacua, on 
which (A^) = and (A^) ~ 0{m^/2)i respectively. We assume that our universe is at the 
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latter, which can be the absolute minimum of the scalar potential for a proper set of the 
parameters. Then, the Majorana mass term of the fermionic component of A^, i.e. m^N'^ 
is generated in the superpotential: 

(N) rriN ^ 0{my2)- (2-3) 

Since we regard the fermionic component of as the extra DM component explaining the 
ATIC/PPB-BETS' observation, we take rriN = "t-dm ~ 2 TeV. 

The non-vanishing VEV {N) breaks U(1)r to the discrete Z2 symmetry, because the 
unit R charge is 1/3 in this model. Since the superfields carrying R = 1/3, 1, 5/3 (0, 2/3, 2) 
become odd (even) under Z2, the remaining Z2 symmetry is exactly identified with the 
R (or matter) parity. In fact, the U(1)r breaking source is the SUSY breaking source 
{F) ~ m^i2Mp ~ (lO^'' GeV)^, which is the VEV of the F-component of a hidden sector 
superfield, and generates the SUSY breaking soft terms in the visible sector. Since the 
R parity of is even, A^ can not be the Majorana neutrino participating in the seesaw 
mechanism. Since the coupling of A^ to the MSSM Higgs doublets is possible, at best, only 
at the high order superpotential, {{N)'^ /Mp)Nhuh(i, it can not also be the extra singlet 
appearing in the "next-to-minimal supersymmetric standard model (NMSSM)" ||l6|| . 

With the terms in the superpotential Eq. (|2.1|), the DM A^ can decays to x^~^^~ via a 
dimension 6 operator induced by a one loop diagram, if moM ^ mgc: 

N — yx + e' + e+. (2.4) 

See the dominant Feynman diagram in Figure 1, which looks similar to "Penguin diagram" 
appearing in the K and B meson decays. As seen in Figure 1, the effective dimensionless 
coupling of e^*Ne^ in the Lagrangian is induced by the loop. It is estimated as 

m3/2{N) 4 

where we set Mj, = Mjr = = Mq = M^. 0{y^) denotes the contributions of the 
dimensionless Yukawa coupling constants, which are assumed to be of order unity at the 
GUT scale. Since the superheavy fields are involved in the relevant Yukawa terms in 



Eq. (2.1), the couplings of the terms ^^NEX" and "XOe*^" do not much evolve with 
energy after the superheavy fields decoupled. [The order of magnitudes of the A^^ coupling 
at the GUT and lower energies are the same because of the small beta function coefficient.] 
Thus, the low energy effective coupling, i.e. Eq. ( |2.5| ), which is obtained by integrating out 
the superheavy particles, is extremely small [< 0(7713^2/-^*)]- If ^^'^ O are in large 

dimensional representations under the other (visible or hidden) non-abelian (gauge) groups 
G, the dimension "AT" can be crucial in Eq. (^). The decay rate of A'' ^ % -|- e~ -|- e+ is 
estimated as 

g'm3/2{N) 
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xO(y«)xA/■^ (2.6) 



where Fat ~ IQ-^^ sec.-^ for moM ~ 2 TeV [> 10 x 0{m^)], ~ lO^^ GeV, 0{y^) ~ 1, 
and J\f = 1. Note that if the dimensionless Yukawa couplings in Eq. (|2.1|) are about 3, 
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Figure 1: Penguin- type one loop decay diagram of N: It is the dominant diagram of N 
X + + . The dimensionless Yukawa couplings are of order unity. 



M* can be slightly heavier upto 10^^ GeV, yielding the same decay rate. The other non- 
abelian (global or gauge) symmetry Q, under which X, and O are charged, would be 
useful in raising M^, higher. For instance, if ^ = flipped SU(5) in the visible sector and 
the superheavy fields are of the SU(5) tensor representation, 7^ = 10 [or ^ = SO(IO) 
in the hidden sector and the superheavy fields are of the SO(IO) vector representations, 
TZ = TZ* = 10], then the circulating fields on the loop are 10 times more (i.e. Af = 10) and 
so the decay rate is 100 times enhanced, compared to the case of the singlets. 



If the selectron e'^ is relatively hght, tit-dm 



> 



rripc 



in Figure 1, and so the two body decay channel, 
becomes enhanced by O(IOO): 



then can be an on-shell particle 
e~ + e"^ opens. Thus, the decay rate 
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(2.7) 



For Fjv giving 10 sec ^, thus, ~ lO"*^^ GeV is not much affected. 

Note that in this model, (anti-) neutrinos and charged leptons heavier than the electron 
are not produced at all from the DM decay. [The muons eventually decay to the electrons 
and (anti-) neutrinos by the weak interaction.] Hence, this model is completely free from 



the constraints on neutrino flux [17| 
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In this model, we have the two DM components, N and the (bino-hke) LSP x- As 
noted in Ref. [^, even extremely small amount of [0(10"^*^) < (n-Af/'T-x)] produce the 
positron flux needed to account for PAMEL/ATIC data, only if the decay rate is enhanced 
by taking relatively light masses of the exotic mediators [10"^^ GeV < < 10"'^^ GeV]. 
Since the other DM component, x can still support the needed DM density /jdm ~ 10~^ 
GeVcm"^, thus, we have extremely large flexibility for the portion of nj^ /n^. 

We have already a TeV scale mass of N . Thus, N can play the role of the well-known 
weakly interacting massive particle (WIMP) such as the neutralino in the MSSM, e.g. if an 
interaction with some other hidden sector fields H and H'^, W D yhNHH^ is introduced. 
Here yh is a Yukawa coupling constant of order unity and the masses of the scalar partners 
of H and are assumed to be of order the electroweak scale. [Then the annihilation cross 
section of N would be in the needed range for explanation of dark matter ((cr|w|) ~ 10~^^ 
cm^s"^).] N could be in a thermal equilibrium state with H, H'^ by exchanging their 
scalar partners down to a proper decoupling temperature defined with hidden sector fields. 
Departure of N from the interactions could leave the relic energy density of order 10~^ 
GeVcm"'^. Alternatively, could be non-thermally produced by decay of hidden sector 
fields. However, we do not specify a possibility, because we have extremely large flexibility 
of n/vr/n^. 

3. Conclusions 

Along the line of Ref. Q, we proposed another SUSY model with two DM components 
{N,x)- A DM could decay to the SM particles only at loop levels, when the exotics are 
the mediator of the decay process. In this model, the extra DM component decays to 
Xe"'"e~ through a dimension 6 operator induced by a penguin-type one loop diagram. Its 
extremely long life time 10^^ sec. required for explaining the observed positron excess is 
caused by the superheavy masses of exotic states mediating the DM decay. Even with 
extremely small amount of N, the positron excess could be explained. This model is easily 
embedded in flipped SU(5), in which e'^ and remain SU(5) singlets. 
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